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Abstract — In this paper, the trust region algorithm was 
used to identify the parameters of the dynamic model of a 
permanent magnet direct current (PMDC) motor, using the 
MATLAB/Simulink Parameter Estimation tool. The 
objective was to estimate the parameters applying the 
square wave, pseudo-random binary sequence (PRBS) and 
random signals in the motor excitation. The obtained 
models were evaluated in open and closed loop, where a 
speed control project was applied using the entire 
eigenstructure assignment. The error between the 
simulated and real curves of velocity and current were 
evaluated by means of the normalized root mean squared 
error (NRMSE). 

Keywords — Control Systems, Parameters Identification . 

I. INTRODUCTION 

Many industrial processes that require position control, 
variable speed, constant torque, rapid acceleration and 
deceleration make use of direct current (DC) motors [1]. 
However, to control systems, it is often necessary to use 
model-based control techniques, e.g., eigenstructure 
assignment, state estimator, robust control, and optimal 
control. 

Thus, to design controllers and satisfy performance 
criteria, it is fundamental to know the dynamic model of 
the system. For control purposes, it is not necessary to find 
an exact mathematical model, but it is necessary to have a 
model that contemplates its main dynamics. 

One of the great problems in controlling DC motors is the 
lack of information regarding the parameters that 
constitute the mathematical model. Often the motor is 
worn out due to its use, which makes its information out of 
date. Thus, a model with imprecise parameters can reduce 
the efficiency of a control system [2]. 

In this context, the study presents a methodology found in 
[2] to estimate the parameters of the dynamic model of a 
permanent magnet direct current (PMDC) motor, using the 
Parameter Estimation tool of the MATLAB/Simulink 
software. For the practical realization of the study, the 
Maxon Motor F2140 shown in Fig. 1 was used. 



Fig. 1: Maxon Motor E2140. 


The work is organized as follows. Section II presents the 
dynamic model of a DC motor. In Section III the trust 
region method is explained. Section IV shows the 
parameters identification using the trust region algorithm. 
In Section V is presented the tracking control system with 
state feedback and entire eigenstructure assignment. 
Section VI presents the results. Section VII concludes this 
paper. 

II. DC MOTOR DYNAMIC MODEL 


Consider the armature circuit shown in Fig. 2. 



The parameter V a is the supply voltage of the armature, R a 
is the armature resistance, L a is the inductance of the 
armature coil and e b is an induced voltage that opposes the 
supply voltage. It is possible to equate the equivalent 
circuit of the motor in (1), using the Kirchoff laws. 

V a -V Ra -V La -e b = 0 .(1) 

The parameters V Ra and v La can be rewritten as a function 
of the current / a , as presented in (2). 

Va ~ Kla ~ La fja ~ e b = 0.(2) 
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B is the viscous friction. It can be concluded that the sum 
of the motor torques should be zero, resulting in (4). 

T e -T a .-T a -T c = 0.(4) 

where T e is proportional to the armature current, T^, is the 
torque generated by the inertia of the motor when the rotor 
undergoes an acceleration and T ^ is the torque originated 
by the effect of the viscous friction: 

.( 5 ) 

d 


Te = K t I a 


T “' =J dt W 


.( 6 ) 

T w = Ba> .(7) 

By substituting (5), (6) and (7) into (4), the differential 
equation in (8) is obtained. 

r\ 

.( 8 ) 


K t I a -JfoM -Ba)-T c = 0,. 

where K t is a torque constant. 

With the electrical and mechanical characteristics equated 
in (2) and (8), they can be rewritten for the armature current 
in (9) and angular velocity in (10). 
d R n K P V n 

.( 9 ) 


T v a T . v a 

'^t a ~~T~ a ~l~ C0+ T~ 

L a L a L a 


K t 


B 


( 10 ) 


d tvf lj i f 

- 7 -u =T Ia ~ 7 w ~T- 
dt J J J 

Equations (9) and (10) can be represented in the state space 
form presented in (11). 

x = Ax + Bu, and y = Cx + Du, .(11) 

where x is the vector of states and the notation of the point 
indicates the time derivative, u is the vector of inputs, y is 
the vector of measured outputs and A, B, C and D are 
constant matrices for a linear system [3]. Then, (12) and 
(13) represent the state space form of (9) and (10). 


V_Ra 

L a 

Kt 
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Ke 

L a 

B 

"7 J 
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r 1 


0 


0 -- 


/ J 


.( 12 ) 


*R 3 & 1 +G ol K].c») 

Equations (9) and (10) can also be represented in the block 
diagram form. For this, it is necessary to apply the Laplace 
transform in both equations considering the initial 
conditions as zero, resulting in (14) and (15), respectively. 

-K e co(s) + TO 


L a s + R a 
K c I a (s) ~ TO 


40 ) = 

" (s) = ]S + B 

In this way, the block diagram in Fig. 3, represents the 
dynamics of the PMDC motor. 


.(14) 

.(15) 


Fig. 3: Block diagram of the PMDC motor [2]. 

III. TRUST REGION METHOD 

Let a scalar function y = In the vector case, f(x t ) 

depends on a vector 9 of N parameters. Then, it can be said 
that the function f{x{) is parameterized by 6 and can be 
represented as in (16). 

y = /Oi,6>).(16) 

After performing N measurements: 

y( 1) =/(*(!), 0).(17) 

y(2) = /Oi(2), 6 ) .(18) 

y(N)=f(Xi(N),d) .(19) 

The parameters vector 6 can be determined by knowing the 
sets {y(l),y(2), and {x i (l),x i (2), ....x^N)} 

[4]. 

According to [4], parameter estimation methods for 
nonlinear systems can be applied to linear systems, this is 
possible because the class of linear systems is a subset of 
the class of nonlinear systems. 

In this sense, the Trust-Region optimization algorithm is 
used to estimate the PMDC motor parameters, through the 
Parameter Estimation tools from MATLAB/Simulink 

software. 

The purpose of the method is to find a vector 6 that 
minimizes the expression in (20). 

N 

F(e) = ^(eO) 2 ) = e T e = ||e(0)|| 2 .(20) 

7=1 

where F(6) is the sum of squared errors and e(j) is the 
error from the attempt to estimate 6 , presented in (21). 

e(J) = y(j ) - /(*;(/). 0)..(21) 

for j = 1, 2, 

The adjustment of the parameters vector 6 using the trust 
region algorithm is done through an iterative process 
starting from an initial vector 6°. 

Most methods of minimization usually choose the 
direction of the step and then decide on its size. However, 
the trust region method determines an upper limit for the 
size of the step to be given and then find the direction to be 
taken [5]. 

The step size is restricted to each iteration /c, by the radius 
of the trust region h k , centered in 6 k as shown in Fig. 4. 
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The direction of the step d k can be determined by 
minimizing the quadratic approximation q k (.) in the trust 
region with center in 6 k and radius h k . Therefore: 


d k E argmin{q k (d ) : ||d|| < h k ], .(22) 

similarly: 

1 

min q k (d ) = F(6 k ) + g T d +-d T Hd, .(23) 

subjected to: 

lldll <h\ .(24) 


where g = F'(0 k ) is the gradient and H = F"(0 k ) is the 
hessian matrix. Thus, after determining d k it is necessary 
to evaluate F(6 k + d k ) in order to verify if 6 k + d k is 
satisfactory [6]: 

F(6 k + d k ) < F(0 k ) .(25) 

In order for (25) to be true, the relation in (26) can be 
assumed. 

6 k+1 = 6 k + d k .(26) 

IV. PARAMETERS IDENTIFICATION USING 
THE TRUST REGION METHOD 

This approach consists of determining the parameters L a , 
R a , K e ,J, B and K t , present in the block diagram of Fig. 3, 
using the MATLAB/Simulink Parameter Estimation tool, 
which is widely used for estimation and optimization of 
model parameters based on experimental data. 

To use the Parameter Estimation tool, it is necessary to 
create a block diagram according to Fig. 3 to represent the 
dynamic model of the PMDC motor. The diagram created 
in MATLAB/Simulink is shown in Fig. 5, and it can be 
seen that K t = K e was adopted to reduce the number of 
variables and simplify the estimation. 



Fig. 5: Block diagram created using MATLAB/Simulink. 


With the diagram created, it is necessary to select in the top 
menu of MATLAB/Simulink the option Analysis and then 
Parameter Estimation to open the tool shown in Fig. 6. Its 
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main window can change depending on the version of 
MATLAB/Simulink, in this project was used the version 
R2015a. 



Fig. 6: Main window of the Parameter Estimation tool. 

As shown in Fig. 6 it is necessary to perform four steps, 
which summarize the procedure of configuration and 
execution of the estimation process: 

1. It is necessary to select the parameters to be estimated 
and optionally the initial conditions and limits. The 
configuration window is shown in Fig. 7 and Table 1 
shows the settings used for each parameter. The value 
of the scale setting can be used to normalize the data. 



Fig. 7: Parameters configuration window. 


Table 1: Initial parameters settings. 


Parameter 

Initial 

Value 

Maximum 

Minimum 

Scale 

B ( Nm. s 
/rad) 

0.01 

1 

0 

1 

] ( Kg.m 2 ) 

0.01 

1 

0 

1 

K e (V.s 

/rad) 

0.01 

1 

0 

1 

La 00 

0.01 

1 

0 

1 

R a 01) 

0.1 

20 

0 

1 
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2. It is necessary to import the experimental data obtained 
from the real system. The armature current I a (t) and 
the speed of the rotor co(t) are the output data. The 
armature voltage V a (t) is the input data. Fig. 8 shows 
that for each of these data it is necessary to provide an 
array of two columns, the first column is a time vector, 
and the second column is the data itself. 



[TimE CurrEnt] 
■gstirnaeao.-lntearatorl :1 fVeloddade'; 
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Initial States 

0ptio nally define initial states f o rthis experiment. 

There are currently no initial states defined for this experiment. 

^selectintia l states 

Parameters 

Optionally define parameters for this experiment. 

Thsra are currently no parameters defined forth is eKperiment. 
a select Parameters 

|~5 Plot a simulate Qrn ^ ok © Help 


Fig. 8: Experiment configuration window. 

3. This step consists in the configuration of the estimation 
options, such as tolerances, methods and algorithms. 
Fig. 9 shows the estimation configuration window and 
Tables 2 and 3 show the settings used in this work. The 
Parameter Tolerance option causes the estimation 
process to end when the parameter variation is less than 
the specified value. The Function Tolerance option 
causes the estimation process to end when the cost 
function variation is less than the specified value. 
Finally, the Maximum Iterations option defines the 
maximum number of iterations allowed in the 
estimation process. 


i lion Options 


General Options Optimijati pn Options Parallel Options 
Optimization method 

Method: Nonlinear least squares_J v Algorithm; Trust-Region-Reflective 


Optimization options 
Parameter tolerance: D.001 
Maximum iterations: 100 

. Use robust cost 

Display level: Iteration 

Restarts; 0 


Funtbcn tbleranct: D.G01 


Cantef Hefp 


Fig. 9: Estimation configuration window. 


Table 2: Optimization method settings. 


Method 

Algorithm 

Cost 

Function 

Nonlinear 

Trust 

Sum of 

Least 

Region 

Squared 

Squares 

Reflective 

Errors 

Table 3: Optimization options settings. 

Parameter 

Function 

Maximum 

Tolerance 

Tolerance 

Iterations 

le-09 

le-09 

100 


After performing the previous steps, the estimation 
process is started. The Parameter Estimation tool 
calculates the cost function by comparing the 
experimental data obtained from the real system with 
data obtained using the model created in 
MATLAB/Simulink. Fig. 10 shows a complete 
estimation and the values of the parameters obtained 
for a given experiment. 



Fig. 10: Complete estimation process. 

This procedure was performed for three different 
approaches, each using an excitation signal. As shown in 
Figs. 11a, lib and 11c, respectively, the used excitation 
signals were square wave, pseudo-random binary sequence 
(PRBS) and random. 

20 -1-1-1-1-.-1-1-1-1- 
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Time (s) 
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Time (s) 


-1 - 1 - 1 — 




0 1 2 3 4 5 6 

Time (s) 


9 10 


Fig. 11: (a) Square wave signal, (b) PRBS signal. 

(c) Random signal. 

V. TRACKING SYSTEM WITH STATE 
FEEDBACK 

The dynamic model obtained for the PMDC motor was 
validated in closed loop through a speed control using a 
tracking system with state feedback, which was also 
applied in the real plant to compare the results. 
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Consider the block diagram in Fig. 12. The open-loop 
system is given by the state space form in (27). 

x = Ax + Bu, and y = Cx = [JT] x .(27) 


The vector w = Ex represents the outputs that are required 
to follow the input vector r that consists of piecewise- 
constant command inputs [7]. 
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V(A + BK) = {v v v 2 ,..., v n ], .( 38 ) 

which are selected in order to obtain the desired time 
response characteristics. The eigenvalues and eigenvectors 
can be related by (39) that can be rewritten as (40). 

[A + BK]Vi = A i v i , .( 39 ) 

[A-A t I B][p.] = 0.(40) 

for i = 1,2 ...,n, where v t is the eigenvector and g t is 
given by (41). 

gt = Kv t .(41) 

To satisfy (40), the vector [vf gJ] T must belong to the 
kernel in (42). 

S(A i ) = [A-A i I B], .(42) 

for i = 1,2 


Fig. 12: Block diagram for the tracking control system 

[71 

This control approach uses an integrator that makes the 
error between the reference and the controlled state go to 
zero when the system is in steady state [7]. 

Therefore, to determine the gains K t and K 2 that allocate 
the eigenvalues of the closed-loop system in the left half¬ 
plane, it is necessary to add in equation (27) the integrator 


dynamics given by (28). 

z = r- w = r- Ex .(28) 

The state space form of (27) and (28) is given in (29) and 
(30). 

GH4 acHM-.^ 

y = [c o][*].(30) 

Where u is given in (31). 

u = K 1 x + K 2 z=[K 1 K 2 \ [*].(31) 


The closed-loop system is given in (32), by substituting 
(31) into (29). 

EMC, .p 2 ) 

Where x' and x' are given in (33). 

x' = (A + BK)x' + B'r, and x' = [*].(33) 


Therefore, the augmented matrices are given in (34), (35) 
and (36). 



K=[K i K 2 \ .(36) 

Then, it is necessary to determine the state feedback gain 
K that stabilizes the augmented system ( A, B ). However, 
the application of the control law is possible if and only if 
the matrices pair ( A, B ) is controllable [7]. 

5.1 Entire Eigenstructure Assignment 
This methodology consists in obtaining the matrix K from 
the selection of the eigenvalues to be assigned to the matrix 
of the closed-loop plant of order n in (33): 

o(A + BK) = {A 1 ,A . . A n ] .(37) 

and an associated set of eigenvectors: 


The notation kerS(Ai) is used to define the space called 
null that contains the eigenvectors [vf gYY in order for 
(40) to be satisfied [7]. Equation (41) can be used to form 
the matrix equality in (43). 

1 0i 9i 9n\ = [Kv 1 Kv 2 ••• Kv n ] .(43) 

The matrix K is obtained through (43) and is presented in 
(44). 

K=[g l g 2 9n][Vl V 2 ••• 

= QV~ 1 (44) 

The eigenvalues can have repeated numbers equal to the 
system inputs. This is because the null space has a 
dimension equal to the number of inputs. In this way, a 
repeated eigenvalue is associated to a vector of the basis of 
the null space. Thus, all the columns of the matrix V remain 
linearly independent and, therefore, the matrix V~ x exists 

m. 

VI. RESULTS 

The parameter estimation procedure was performed for 
each excitation signal: square wave, PRBS and random. 
The curves resulting from the estimates are shown in Figs. 
13, 14 and 15, respectively. 



1 i — i i-—r-- 't- 



0123456TB9 10 


Tims (s) 

Fig. 13: (a) Real and simulated speed for the square wave 
signal, (b) Real and simulated current for the square 
wave signal. 
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Fig. 14: (a) Real and simulated speed for the PRBS 
signal, (b) Real and simulated current for the PRBS 
signal. 
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Fig. 15: (a) Real and simulated speed for the random 
signal, (b) Real and simulated current for the random 
signal. 


By analyzing the results, it can be verified that the 
parameters identified for the PMDC motor are very close 
to the actual parameters provided by the manufacturer, as 
shown in Tables 4 and 5. 


Table 4: Estimated parameters. 



Estimated Parameters 

Parameter 

Square 

Wave 

Signal 

PRBS 

Random 


Signal 

Signal 

B ( Nm.s/rad ) 

2.33e-06 

1.73e-06 

1.75e-0.6 

J ( Kg.m 2 ) 

2.32e-06 

2.42e-06 

2.39e-06 

K e C V.s/rad) 

0.027439 

0.028006 

0.027971 

La 00 

0.0015168 

0.001261 

0.014311 

Ra 01) 

10.223 

10.489 

10.639 


Table 5: Parameters provided by the manufacturer. 


Parameter 

Manufacturer’s 

Parameters 

B ( Nm.s/rad ) 

not provided 

J (Kg.m 2 ) 

2.30e-06 

K e ( V.s/rad ) 

0.0278 

L a 00 

0.00127 

Ra (LI) 

10.7 


To evaluate the quality of the estimated parameters, three 
validation signals were applied in the real motor and in the 
simulated models as shown in Figs. 16, 17 and 18, 
respectively, the step, sinusoidal and triangular signals. 
Table 6 shows the results of the comparisons between the 
curves obtained experimentally and through the 
simulation. 

In Table 6, a higher value of the normalized root mean 
square error (NRMSE) criterion indicates a greater 
proximity between the compared curves and, therefore, 
indicates a better result. Thus, an NRMSE value of 100% 
indicates that the curves are the same. 


■Real 
- Simulated 


0 1 2 3 4 5 6 

Time (s) 


8 9 10 


0) 

O -0.5 - 





"Real 

- Simulated 



2 3 4 5 6 

Time (s) 


8 9 10 


Fig. 16: Validation response by means of a step input 
signal. 




Fig. 17: Validation response by means of a sinusoidal 
input signal. 
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Fig. 18: Validation response by means of a triangular 
input signal. 


In Table 6, the NRMSE for the armature current l a (t) 
presents small values, this is because there were 
nonlinearities in the model that were not considered, e.g., 
the Coulomb’s friction and the dead zone of the motor (i.e. 
for a range of applied voltages, the motor does not rotate). 
These nonlinearities are significant when the motor passes 
through the zero speed point when rotating in both 
directions [2]. However, as this work assumes that the 
mathematical model of the PMDC motor is linear, the 
results are considered satisfactory. In addition, it can be 
verified that the results of the models identified for the 
square wave, PRBS and random signals were close and it 
is not possible to indicate the best result. 


Table 6 : Comparison between the results obtained 
experimentally and through the simulation. 





Excitation Signals 

NRMSE (%) 

Used for Estimation 

Square 

Wave 

PRBS 

Rando 

m 


Step 

<w(t) 

96.804 

5 

97.09 

17 

97.015 

2 

in 

13 

I ait) 

63.232 

8 

59.74 

82 

59.121 

9 

= 

.£f 

35 

= 

Sinusoid 

(O(t) 

98.256 

0 

98.14 

91 

98.151 

0 

j © 

’-£2 

a 

2 

al 

/«(t) 

83.633 

5 

82.79 

29 

82.904 

6 

a 

> 


<w(t) 

97.419 

97.77 

97.775 


Triangul 

0 

56 

5 


ar 

I ait) 

48.907 

2 

45.17 

23 

45.508 

4 

Mean (%) 

81.375 

5 

80.12 

16 

80.079 

4 


square-wave excitation signal was chosen. Then, the 
control method of entire eigenstructure assignment was 
applied to the PMDC motor. The eigenvalues A lf A 2 and A 3 
were chosen empirically in order to ensure that the time 
response presented the desired characteristics and a 
satisfactory result, without presenting overshoot. In (45), 
(46) and (47) are presented, respectively, the eigenvalues 
A ± , A 2 and A 3 . 


A ± = -22.4471.(45) 

A 2 = -260.5043.(46) 

A 3 = -320.2260.(47) 


From the eigenvalues in (45), (46) and (47) the gains K x 
and K 2 presented, respectively, in (48) and (49) were 
obtained. Where K 1 represents the gains of the states and 
K 2 represents the tracking system integrator. 


= [9.3096 0.0151].(48) 

K 2 = [0.2401].(49) 


After the gains K 1 and K 2 were determined, three 
reference signals were applied for validation of the control 
and comparison of the curves obtained experimentally and 
through the simulation. Figs. 19, 20 and 21 shows the 
control responses, respectively, for the step, square and 
sine wave validation signals. 

600 


-o 400 - 


3 0 

o 


-0.5 




■a;(t) Reai 
•u;(t) Simulated' 


-Reference 


0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Time (s) 



3 4 5 


Time (s) 


Time (s) 


Fig. 19: Control response for the step signal. 



Time (s) 




Fig. 20: Control response for the square wave signal. 
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In this way, only one of the obtained models was used for 
the closed-loop evaluation, the model obtained with the 
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!- W F pffl 

- Ia(t) Simulated 

-0.4 I——I .20 -----*- 


10 15 20 

Time (s) 


10 15 20 

Time (s) 


Fig. 21: Control response for the sinusoidal signal. 


The Table 7 presents, for the control responses, the 
comparison between the curves obtained experimentally 
and the curves obtained through the simulation. It can be 
verified that for the control project, the simulated model 
satisfactorily described the actual PMDC motor, since the 
speed responses for all references were close to the 
experimental responses. 

Table 7: Comparison between the real and simulated 
curves of the control responses. 


Reference Signals 

NRMSE (%) 


6>(t) 

93.2833 

Step 

Ia(.t ) 

26.2284 


u(t) 

86.9883 

Square 

Wave 

6>(t) 

94.7700 

/a(t) 

41.2453 

u(t) 

91.4634 



98.2582 

Sinusoidal 

Ia(t) 

43.4177 


u(t) 

96.3833 

Mean (%) 

74.6708 


VII. CONCLUSIONS 

In this work, the trust region algorithm was applied to 
determine the parameters of a PMDC motor. Three 
excitation signals were used and the efficiency of the 
method was verified for the three cases, all the identified 
models satisfactorily represented the behavior of the 
PMDC motor. 

In order to arrive at the results the models were validated 
in open-loop applying the step, sinusoidal and triangular 
signals, and comparing the curves of velocity and armature 
current, obtained experimentally with the curves obtained 
in the simulation. 

Then, these models were evaluated in closed-loop, where 
a tracking system with state feedback was applied for 


speed control. The experimental responses were 
satisfactorily close to the simulated responses. 

It is concluded that, in practice, the estimation of PMDC 
motor parameters by the trust region method with the 
Parameter Estimation tool is simple to use and efficient. 
The next step will be to compare this identification method 
with traditional methodologies and also using Kalman 
Filter Unscented. 
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